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Abstract 

For the special case of compact QED in (2+1) dimensions, we calculate the non- 
Gaussian vacuum wave-functional to second order in the monopole fugacity and 
obtain the effective photon mass. Our method presents some hope for understanding 
the connection between variational and systematic approaches to understanding the 
non-perturbative wave- functional. 



1 Introduction 



In recent years there has been a revival of interest in Hamiltonian methods for Yang-Mills 
gauge theories. These studies have for the most part fallen into two groups: one group 
tries to solve the functional Schrodinger equation in some sort of systematic expansion 
[l]-[9], the other makes use of a Gaussian ansatz for the wave-functional, and a variational 
principle that minimises the vacuum energy on the Gaussian state [10]-[14]. Both of these 
approaches have led to interesting insights into the vacuum structure of non-abelian gauge 
theories, but in both cases many important questions remain unanswered. 

Systematic methods run into a number of technical difficulties that have been resolved 
in some cases but not in general. One of the first problems is to rewrite the Hamiltonian 
in terms of suitable gauge-invariant variables. This must be done in such a way that the 
relevant systematic expansion is well-defined and renormalisable. Certain expectations, 
such as the existence of a derivative expansion of the wave-functional at large distances 
[15] are frustrated by complications arising from the existence of topologically non-trivial 
classical field configurations [7, 16]; some choices of gauge-invariant variables can take 
them into account implicitly, or one can perform an explicit expansion around classical 
solutions with non-trivial toplogy. Dynamically generated mass parameters emerge from 
the analysis in various ways [7, 11]. 

The variational approach sidesteps many of these issues, at the cost of restricting to 
a (gauge-projected) Gaussian state. This allows many concrete calculations to be per- 
formed; on the other hand, though it is reasonable to expect that the Gaussian state 
captures the important physics, this is difficult to justify. It would be desirable to "ele- 
vate" the results of the variational method to the full non-Gaussian theory, or at least to 
understand the relation between the two. 

In this paper we will achieve this for the simple case of compact QED in (2-1-1) di- 
mensions, extending the analysis of [11] to the full non-Gaussian theory. This model 
exhibits many interesting features such as dynamical mass generation and confinement, 
and though the mechanisms for these features are special to (2+1) dimensions, many 
general features of our method are more widely applicable. Our analysis makes use of a 
systematic expansion in the monopole fugacity, which is closely related to the dynami- 
cally generated photon mass. Although the dynamically generated mass is of the order of 
the monopole fugacity, to calculate it in the non-Gaussian theory requires an analysis to 
second order in this parameter. Results to higher order are obtainable in principle. 

The variational principle is essentially trivial, and the meat of the calculation is in 
the effective monopole dynamics, but the former nevertheless allows us to determine 
the contributions of higher n-point functions to the propagator and solve the functional 
Schrodinger equation order by order. Normalisability conditions on the wave-functional 
also play a role. 

The ultimate aim of our investigations is to make progress towards Hamiltonian meth- 
ods for the study of QCD and Yang- Mills in (3+1) dimensions. It is our behef that though 
such methods are currently rudimentary, they can be developed to the point where such 
methods can rival the explicit calculations available, for example, in lattice gauge theory. 
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2 Free QED 



To begin with consider free non-compact electrodynamics in (2+1) dimensions. The 
Hamiltonian is 

H=\{B^ + E^l (1) 
where B = eijdiAj and the electric field operator is represented as 

B, ~ A. (2) 

If we make a Gaussian ansatz for the vacuum wave-functional 

^B] = exp{-lBCB}, (3) 

(here we have made use of a matrix notation, so that BCB = J (fxd'^yB{x)C{x — y)B{y), 
etc.) then the Schrodinger equation iJ^ = E'^ implies that (in momentum space) 

-p'Cipf + 1 = 0, ElV=\j ^yC{p). (4) 

Here E is the vacuum energy and V is the volume of space. From (4) we immediately 
obtain C{p) = 1/p. 

We can also derive this result by means of a variational principle: according to the 
standard prescription for calculating expectation values we have 

^p2^ I DAB^^B]^*[B] _ 1 /• d^p 1 

while 



J DAi-^[B]-^*[B] 2 J 47r2' 

Minimising the expectation value of the Hamiltonian then gives C{p) = 1/p as before 
(and the virial theorem for harmonic oscillators is satisfied since (Ef) = (B"^)). 

There is a third way of deriving this result that does not depend on making a Gaussian 
ansatz. If we assume a more general form for the vacuum wave-functionaP 

^[B] = exp{-lBCB - ^C,B' - ^C^B^ -...}, (7) 
then the Schrodinger equation becomes 



= _c(p) V + iC4(p)p2 + 1 

; (8) 



^In matrix notation 64-8^ = / (Pxi...(PxiCi{x\,X2,Xz,Xi)B{x^)B{x2)B{x^)B{x^), etc. 
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In (8) C4(p) is the Fourier transform of C4(a;) = / cPyCA{Q, y, y, x), and we have assumed 
translational invariance. The expression for the vacuum energy density is unchanged, 
while the order part of the Schrodingcr equation acquires a contribution from the 
four-point function. The vertical dots indicate the existence of further equations at order 
5^ and so on that will not concern us much for the moment. 

In addition to the equations (8), the vacuum wavefunctional must satisfy normalis- 
ability conditions C > 0, C4 > 0, etc. If we minimise the vacuum energy density subject 
to these normalisability conditions we get C(p) = 1/p as before, while C4 and all other 
n-point functions vanish. Note that the energy minimisation determines the contribution 
of the 4-point function to the 2-point function, allowing us to solve (8) order by order in 
B^. This is a general feature of our methods. 

3 Compact QED 

The compact theory differs from the non-compact one by the presence of Dirac monopoles, 
which play the role of instantons in three space-time dimensions. As usual, we can repre- 
sent the vacuum as a sum over instanton sectors 

nB] = Y.'^n[Bi (9) 

n 

where represents the vacuum in the n-monopole sector. Equivalently, we can require 
the vacuum to be invariant under the action of the vortex creation operator 

nx)=exp{i/A^|[^iJ,w}, (10) 

so that the vacuum wave-functional is a superposition of wave-functionals of the non- 
compact theory with an arbitrary number of vortices and anti-vortices at any spatial 
point [11]: 

= n / ^'^aAV(xJ n / d-'xpK'V^xpmB]. (11) 

^+-''^-- n+,n_=Oa=l 0=1'' 

We have introduced an ultraviolet cutoff A. If the theory is regarded as the unbroken 
sector of the Georgi-Glashow model, then A is related to the charged vector boson mass 
Mw', the details depend on the UV dynamics [17]. The reason for the aforementioned 
equivalence is as follows: inner-products of the form (5) are regularised by point-splitting, 
and monopoles that live on the quantisation surface have Dirac lines that can intersect just 
one of the copies of the wave-functional in the inner-product, inducing a vortex, or neither. 
Thus the vacuum, summed over monopole sectors, is invariant under multiplication by 
the vortex operator. 

The vortex/monopole is represented on the wave-functional by a shift of y(5^(a;) in 
the magnetic field B{x). Introducing the monopole density p{x) = — Yl,a,/3{^'^{^ — Xa) — 
5^{x — Xj})), the sum over monopoles in the non-Gaussian wave-functional (7) gives 
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MB] = Y: exp{-i(S + p)C{B + p)- iC,{B + pf - ^C,{B + pf - . . .}. (12) 
p 

We have normalised the wave- functional to ^[0] = 1 with the factor Z — Z^^^fp]. The 
vacuum energy density is then 



E/v = \z-' E / [p'c{p) + [\p'p'cM - p\pc{p)f}\ + oip')] (13) 

where pC{x) — / (Pyp{y)C{x — y), etc. 

Expectation values of the monopole density are often calculated in the effective low 
energy theory, which takes the form of a sine-Gordon theory [11, 17]. They can also be 
calculated directly by expanding the wave-functional in monopole sectors; thus 

Z = Y.m[p\ = l + 2Vz + ... (14) 
p 

where z — A^exp{— ^C(O) — |^C4(0, 0, 0, 0) — . . .} is the monopole fugacity, and we 
have inchidcd only the dominant one (anti-)monopole sector explicitly. Note that the 
expression for the monopole fugacity is to be understood in a UV regularised sense, so 
that C(0) ~ C(l/A), etc. It can be shown that all volume factors cancel in the expression 
(12) (logarithmic divergences appear in the wave-functional at higher orders, and can be 
removed by normal-ordering). This kind of structure is familiar in mass perturbation 
theory [18, 19]. 

To first order in z we have 

(p(pi)...p(p„)) =Z-i^p(pi)...p(p„)*[p] = — ^52(pi + ...+p„). (15) 

p 3 

Now to this order the expression (13) for the energy density can be simplified by using 
the equations (8). We do not need to consider 0{z) corrections to these equations, as 
they will lead to 0(2;^) corrections to the vacuum energy. The order p^ terms in (13) give 
an order z term, and the order p^ and higher terms cancel, giving 



P C[p) 

9^ 



p ■> "t "" L y J 

This is the same result we would have obtained if we had assumed that the wave-functional 
tends smoothly to the wave-functional of the non-compact theory as 2; — > 0, so that at 
order 1, C{p) — 1/p and C4(p) = 0, etc. 

To proceed further it is useful to expand the logarithm of the wave- functional in powers 
of the magnetic field. So we write 



In* = In ( e-i^^^ ~ ^^4^^ ~ ■ ■ ■ 
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±^ J <fxC{xo-x)B{x)+^ J cPx(PyC4(xo,xo,x,y)B(x)B{y)+.. 



X [1 + Z I d Xo e " ^ " ' ' ' 9 



+ ... 



= -lBC'^ffB + jfCl^^B'' + 0{B') + 0{z'). (17) 

Here the sum over ± corresponds to the one-monopole/one anti-monopole sectors and 
we have exactly the same simphfication that we observed in the vacuum energy, giving a 
shift C{p) — >■ C^^^ip) = C{p) — ^-4- For higher orders in B we find 



(Pi, . . . ,P2n) - C,^ {2n-2)lg-n- ^ 

Writing the equations (8) in terms of C^^^, C^/^^ ^ etc. we have (to first order in the 
monopole fugacity) the same equations and normahsabihty conditions as before, and 
hence the same solution C^-'^-^(p) = = for n > 1. Thus C(p) = ^ + etc. 

Our failure to see mass generation at this order does not come as so much of a surprise 
if we re-examine the variational Gaussian calculation of [11]. There the order z ~ rr? 
contributions from the electric and magnetic energy densities cancel against one another. 
To see the mass generation in our formalism, we will have to extend our analysis to order 
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Z 



4 The second-order calculation 

If we include two (anti-)monopole configurations in the calculations of the last section, 
we find the two monopole contribution 

(p(x)p(y))2 = z^ (±e±^(^-^) + ?>\x - y) d^^w (e^^^"') - l)) . (19) 

Here the sign is ± as the monopoles have equal or opposite charge, and in the second 
term we have a partial cancellation of the infrared divergence due to the normalisation 
factor Z, though a logarithmic divergence remains. The monopole density in the vacuum 
is low for z <^ 1 so that (11) is the partition function for a dilute gas, in which monopoles 

are assumed to be widely separated, hence the cutoff in the integral. 

Putting this result into (17) we find the order z^ contribution to C^-^^{x — y): 

d^Xid^X2-. :-. r siuh 




\xi-x2\>L \x-xi\\X2~y\ Vfi'^|a;i - a;2| 




(cosh (^j-lj I . (20) 
In the dilute gas approximation we can take sinh ( -jr^^ — r) ~ -^r-^ — r and evaluate the 

° -^-^ VS Fl-^2|y g''\x\-X2\ 

first integral in (20) as This corresponds to a mass term in the two-point function; 

the puzzle is that the mass w? ~ ^ appears to be of order z"^, whereas we know from [17] 
that the effective photon mass is of order z. 
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Again, the problem is that the topological expansion makes it difficult to read the 
physical mass off from the wave-functional directly. What we need to determine the 
physical mass is a vacuum expectation value like {B{x)B{y)) , which can be read off 
from the functional Fourier transform of the non-Gaussian wave-functional that we have 
constructed. 

To obtain the functional Fourier transform we treat the four-point and higher functions 
as perturbations of the quadratic wave-functional. The Schrodinger equation again has 
the solution C^^^{p) — 1/p, C^l/ = for n > 1, and we can write the wave- functional in 
terms of the field Ai (with vector indices absorbed into the matrix notation) as 

■^[A,] = eM-\ACA - \C^A^ - ^^6^6 -...}, (21) 

where for example Ciji) = Vp'^C{p) with V a transverse projector. For our purposes we 
can now forget about the topological expansion, since VEV's of the magnetic field operator 
are unaffected by the presence of monopoles [11]. We are interested in the quadratic part 
of the Fourier transformed wave-functional, from which we can read off the physical mass. 

The latter receives contributions from tadpole diagrams like / (Pxi . . . (Px2nC2n{xi, ■ ■ ■ , X2n) 
C~^{x — Xi)C~^{x2 — Xs) . . . C~^{x2n — y), etc. which, when summed, give the final result 
for the photon mass 

'"^~7''"'(-v) = 7^- '''' 

This result agrees with [11, 17], and is of first order in the monopole fugacity, as expected. 
Interestingly, the dynamical mass generation is given here by the condensation in the 
vacuum of a monopole/anti-monopole pair. 



5 Conclusions 

For the case of compact QED in (2+1) dimensions we have shown how the mass genera- 
tion first demonstrated by Polyakov [17] and reproduced in the Hamiltonian formalism in 
[11] by means of a gauge-projected Gaussian ansatz, can be seen in the full non-Gaussian 
wave-functional. There is a variational principle implicit in any solution of the Schrodinger 
equation, and our method provides a link between variational methods based on a Gaus- 
sian ansatz, and attempts to solve the functional Schrodinger equation exactly in some 
non-pcrturbative expansion. 

In our analysis the variational principle and normalisability conditions on the wave- 
functional played a role, but the mass generation was seen to be a feature of the effective 
monopole dynamics. In general, explicit expansion around topological solutions may be 
the best way to tackle a systematic Hamiltonian analysis of non-abelian gauge theories. 

The extension of these functional Schrodinger methods to fermions and supcrfields 
is straightforward. The study of supersymmetric theories along these lines is strongly 
suggested by the non-renormalisation theorems afforded by such theories; these simplify 
the quantisation considerably, as well as facilitating the interpolation between different 
perturbative regimes. 
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